ON LERAY'S PROBLEM FOR ALMOST PERIODIC FLOWS 
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Abstract. We prove existence and uniqueness for fully-developed (Poiseuille-type) flows in semi- 
infinite cylinders, in the setting of (time) almost-periodic functions. In the case of Stepanov almost- 
periodic functions the proof is based on a detailed variational analysis of a linear "inverse" prob- 
lem, while in the Besicovitch setting the proof follows by a precise analysis in wave-numbers. 

Next, we use our results to construct a unique almost periodic solution to the so called "Leray's 
problem" concerning 3D fluid motion in two semi-infinite cylinders connected by a bounded reser- 
voir In the case of Stepanov functions we need a natural restriction on the size of the flux, while 
for Besicovitch solutions certain limitations on the generalized Fourier coefficients are requested. 
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1. Introduction 

We consider the motion of a viscous fluid in semi-infinite cylindrical pipes, with an assigned 
(time) almost-periodic flux. The results are aimed to find solutions of the so-called "Leray's 
problem." Moreover, this work can be considered as a intermediate step towards the analy- 
sis of (deterministic) statistical solutions for the flow of Poiseuille-type, which is the object of 
our current and ongoing research. We recall that Leray's problem (which seems to have been 
proposed by Leray to Ladyzhenskaya [ ]) is that of determining a motion in a region with 
cylindrical exits, subject to a given flux, and tending to the Poiseuille solution in each exit. More 
precisely, let be given a connected open domain O c R'^ made of a "reservoir", a bounded and 
smooth open set Oq, with two cylindrical exits Oi and O2. These two semi-infinite exits (pipes) 
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are described in coordinate systems directed along the axis as 

Oi = Di X R+, 



where the smooth cross sections D^, i = 1, 2, may be possibly of different shape and measure. 
We denote by z e the axial coordinate in both cylinders. Pioneering results in the stationary 
case are those of Ladyzhenskaya [18] and Amick [2]. See also the review in Finn [ ]. The 
extensive literature on the stationary problem is recalled for instance in [ , ] for the linearized 
and full Navier-Stokes problem, respectively. More recently the problem of motion in pipes has 
also been addressed in the time-evolution case, see Ladyzhenskaya and Solonnikov [20] and 
also the review in Solonnikov [ ]. In the last decade Beirao da Veiga [ ] and Pileckas ['^ '] 
gave new contributions to the study of the time-dependent problem with assigned flux, and the 
special role of the pressure has been also emphasized by Galdi and coworkers [14, 13]. 

In [ ] Leray's problem has been considered in the context of time periodic flows, especially 
in view of application to the study of blood flow and we recall that the role of blood flow in 
mathematical research has been put in evidence by Quarteroni [25]. We also stress that the 
(non-trivial) explicit solution introduced by Womersley is periodic, in some sense generalizes 
the Poiseuille flow, and has been discovered in the study of physiological flows. Since the heart 
is pumping with a flux which is not periodic, but a superposition of possibly non-rational fre- 
quencies, this suggests also to study the problem in the setting of almost periodic functions. This 
work has been originated by the inspiring results in [ ] and especially from Remark 3 therein: 
The independence of the various constant on the period of the flux let the author suggest about 
the possible extension to almost periodic solutions. The problem nevertheless requires a precise 
functional setting in order to detect the largest class of almost periodic functions to be employed. 
Moreover, it seems that the very-nice proof based on Fourier series in [ ] cannot be directly ap- 
plied to the new setting and in addition new difficulties in treating the nonlrnearities arise when 
almost periodic functions are employed. This leads us to propose two different approaches in 
two different functional settings. We finally remark that, in addition to early results of Foias [ ], 
the approach via almost-periodic functions finds wide applications in fluid mechanics (see for 
instance the recent paper by Gerard- Varet and Masmoudi [ ]). 

1.1. Setting of the problem. The problem we wish to solve is to find a (time) almost-periodic 
solution of the Navier-Stokes equations 



such that u converges in both pipes as |z| — > oo (in a sense we shall make clear later) to the 
solution of the Poiseuille-type problem. For clarity we recall (cf.[ , § 2]) that by solution of 
the Poiseuille-type problem (of fully-developed flow) we mean a solution of the Navier-Stokes 
equations such that, in a reference frame with z directed along the axis of the pipe and x :— 
(xi, X2) belonging to the orthogonal plane, is of the form 



(1.1) 




xeO, teR, 
X e O, teR, 
on 90, teR, 



u(t,x,z) = (0,0,w(t,x)) 



and 



P(t,X,z) = 7T(t,X,z) +po(t). 
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Here po[t] is an arbitrary function and in addition the flux condition is satisfied 

w(t, x) dx — f (t), 

. D 

for some given function f , where D is the section of the domain. 

The Poiseuille-type ansatz implies that Navier-Stokes equations reduce in the semi-infinite 
pipes Oi, i = 1, 2, to the following equations 

9tw^ - vA^w^ + 9^p^ = 0, (x,z] G Oi, t e R, 

9xiP' = 9x.p' = (x,z) e Oi, te R, 

w^(t,x) =0, X e 90i, t e R, 

where denotes the Laplacian with respect to the variables Xi and X2. From the first equation 
it follows that 9zP^ is independent of z. From the second equation, we also obtain that p^ is 
independent of x, hence p^(t, x, z) = — 7T^(t] z + Po(t). Since the term Ppft) does not affect the 
velocity field, we may assume that the pressure has the form p^(t, z) = — 7t(t)^ z. Moreover, the 
dependence of on the space variables Xi and X2 allows us to consider a problem reduced to 
the cross section Di of Oi, and the flux condition is w^(t, x) dx = f (t). This implies that we 
have to study in each pipe the following problem (called in the sequel the "basic flow"): Find 
(w^(t,x),7T^(t]] such that 

{9twMt,x) - vAwHt,x] = 7t^(t], X e Di, t e R, 

w^(t,x) =0 X e 9Di, t e R, 

Jq. w^(t,x) dx = f(t] teR, 

showing (under suitable assumptions) that if f is almost periodic, then the couple (w, tt] is 
almost-periodic, too. 

We observe that, contrary to the stationary problem where the same approach gives the well- 
known Poiseuille solutions, the solution of the time-dependent motion is more complex for the 
determination of the non-constant pressure (Observe that the classical Poiseuille solution is that 
obtained for circular pipes, but nevertheless in more general domains the same approach gives 
corresponding results). In our problem for the basic flow we have two scalar unknowns and 
two equations, but contrary to the classical problems in fluid mechanics one cannot get rid of the 
quantity tt by means of projection operators. The problem we have to solve can be considered 
as an inverse problem. Moreover the problem cannot be treated with the standard variational 
tools in a direct way. We can write a single equation (the "elimination" of rc^ is obtained by 
taking the mean value over Di) obtaining 

9tw^ - vAw^ +ySu ^■^^ = f X e Di, t e R, 
w^(t,x) =0, X e 9Di, t e R, 

and the latter equation makes easy to understand why some knowledge also of the derivative 
of the flux will be needed in order to solve (2.2). Moreover the usual energy- type estimates 
obtained by testing with w\ or with — Aw^, and with w\ are not-conclusive when applied to this 
problem. In particular, the lack of coercivity prevents from a direct application of the standard 
techniques employed for parabolic problems, see [ , Sec. 3]. This particular issue has been 
addressed in two different ways by Beirao da Veiga [ ] (periodic case) and Pileckas [ ] (given 
smooth flux). 

Even if we generalize to the almost-periodic setting the periodic results obtained in [ ], in 
the first part of the paper we will mainly follow and suitably adapt the approach of [ ]. In 
particular we give special emphasis to the solution of (1.2) since this represents one of the main 
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technical difficulties. The nonlinear problem is then treated by means of perturbation argu- 
ments in a more or less standard way. We want also to point out that in the huge literature on 
almost periodic solutions we find particularly inspiring (for the choice of Stepanov functions as 
suitable for our problem) the paper by Marcati and Valli [ ] concerning compressible fluids. 

In the second part of the paper we consider the problem in the larger class of Besicovitch 
almost periodic solutions with an approach which is more in the spirit of Fourier analysis. We 
give a different proof of the existence of the basic flow which also covers the H^(R) case and 
provides an alternative proof of [ , Thm. 1], when restricted to a time-periodic flux. The fully 
nonlinear case needs, besides the natural assumption of large viscosity, an additional assump- 
tion of regularity on the flux (see (3.12)) which accounts of the technical difficulties of this case, 
due essentially to the non-local (in time) quantities that are used, see Section 3.4.7. 

Remark 1.1. For its variational formulation and the use of energy estimates the problem seems 
to be naturally set in Hilbert spaces and this is not well fitting with the classical continuous 
(Bohr) spaces of almost periodic functions. A suitable choice of the spaces represents then a 
fundamental starting point. We are presenting two different proofs in two different settings, 
since they are substantial different and the assumptions we make on the flux are of very differ- 
ent nature. In the first part we deal with Stepanov a.p. functions, and the setting is much similar 
to the classical variational one for evolution partial differential equations. In the second part we 
deal with Besicovitch a.p. functions and the proof use analysis in wave-numbers. We also point 
out that while the linear problem can be also treated in a unified way, for the nonlinear one the 
differences in the functional setting imply special assumptions on the size of the flux and on the 
Fourier coefficients, respectively. 

Plan of the paper: In Section 2 we consider the problem under the condition of a Stepanov 
almost periodic flux. After recalling the main definition we give a complete solution of Leray's 
problem, with the natural (in space dimension three) restriction of a large viscosity. As by- 
product of our results, we also prove existence in the case of H^(R) fluxes. In Section 3 we 
consider the problem in the larger class of Besicovitch almost periodic solutions and we prove 
existence for the basic flow, together with existence for the nonlinear problem under suitable 
restrictions on the flux. 

2. Leray's problem in the framework of Stepanov a.p. functions 

Here we introduce a functional setting in which it is possible to extend the result of [^] to 
almost periodic solutions. 

2.1. Functional setting. The problem of almost periodic solutions of partial differential equa- 
tions has been studied extensively in the last century, starting with the work of Bohr, Muck- 
enhoupt, Bochner, and Favard and many others. See the review in Amerio and Prouse [ ], 
Besicovitch [ ], Corduneanu [ ], and Levitan and Zhikov [ ]. 

In the sequel we will use the standard Lebesgue and Sobolev spaces H'^ = W*'^. For 
simplicity we also denote by || . || the L^-norm. We will use the symbol C to denote a generic 
constant, possibly different from line to line, depending on the domain and not on the viscosity 
V or on the flux f. Next, given a Banach space (X, || . ||x) we denote by UAP(R; X) the space of 
almost periodic functions in the sense of Bohr-Bochner. We recall that a function f e C°(R; X] 
is almost periodic if and only if the set of its translates is relatively compact in the C"(R; X)- 
topology (observe that if C°(R;X) denotes the space of continuous bounded functions, then 
U AP(R; X) C (R; X)). In the context of weak and strong solutions to partial differential equa- 
tions it is probably better to work with a more general notion of almost-periodicity, given for 
functions f £ LP (R; X), which is suited to deal with distributional solutions. 
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Definition 2.1 (Stepanov p-almost periodicity). We say that the function f : R — > X is Stepanov 
p-almost periodic (denoted by f e SP(R;X)) if f e Lj^^(R;X) and if the set of its translates is 
relatively compact in the Lj|jg^(R; X) topology defined by the norm 



r ft+i 



KIIlp 



(R,X) 



sup 
teR 



When \) — 2 we say simply that the function f : R ^ X is Stepanov almost periodic. 

We will give the main result by using fluxes belonging to this class, together with their first 
derivative. However, in the second part of the paper we will consider also a wider class of 
almost periodic functions: functions almost periodic in the sense of Besicovitch. Further gener- 
alities (not needed in this section) on almost periodic functions are given in Section 3.1. 

A first main result that we will prove concerns the existence of the "basic flow" problem in 
this framework. 

Theorem 2.2. Let be given a smooth, connected, and hounded open set D c R^ and let he given f such 
that f, f e S^(R). Then, there exists a unique solution (w, tt) of (2.2) such that 

Aw,wt e S^IR; L2(D)], 

Vw G S2(R; l2(D)) n CS(R; L^fD)), 

7te S2(R), 



and 



(2.1) 



sup 
teR 



v||Vw(t)|| 



t+i 



(v2||Aw(s)|l2 + ||wt(s)|p + |7t(s)P) ds 



s; c(v + i + 



Remark 2.3. The result concerning the linear problem for the basic flow holds true in any space 
dimension. 

This allows to obtain in a rather standard way the following result for the Navier-Stokes 
equations. 

Theorem 2.4. Let O as in the introduction and let he given f such that f, f e S'^fR]. There exists 
Vq — volf, O) > such that ify > Vo there exists a unique solution u of (1.1) such that 

ue S2(R;H'^(0)) for all s<2 

and u converges to a Poiseuille-type solution in each pipe, as \z\ +oo. 

Remark 2.5. The restriction on the viscosity is not surprising and is common to several results 
concerning the three-dimensional Navier-Stokes equations. This is also observed in { ^-, Ch. XI] 
since the existence of a flux carrier that can be absorbed by the dissipation for any positive 
viscosity is generally not known for cylindrical domains. This imposes (also in the stationary 
case) limitations on the size of the flux, in terms of the viscosity. 

We observe that also in the time-periodic case [ ] largeness (in terms of data of the problem) 
of the viscosity is required. Nevertheless, the results in [ ] concern weak solutions and unique- 
ness is not stated. On the other hand in [ ] there is no restriction on the viscosity, since special 
"two-dimensional-like" solutions are considered. 
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2.2. Construction of the solution of the "basic flow". In this section we give a detailed analysis 
of the existence of an almost periodic basic flow and a complete proof of Theorem 2.2. The 
problem is the following: given f, f G S^(R) find a Stepanov almost periodic solution of 

{9tw — vAw = 7t, X e D, t e R, 

w(t,x) =0, xe 9D, te R, 

Jpw(t,x)dx = f(t) teR, 

Remark 2.6. It is easy to check that one can analyze the slightly more general problem where (2.2) 
is replaced by 

{9tw + yAw = 7te, teR, 
(w(t),e)H =f(t] teR. 

with an imbounded, linear, and with compact inverse operator A on the Hilbert space H with 
domain D(A) and e e H with e ^ D(A) is given. Under suitable assumptions on A, the same 
procedure that we will employ can be used, see also [ ] . The same remark holds also for the 
results of Section 3. 

We start by solving the following initial-boundary value problem in the unknowns (w, 7t), 
'9tw-vAw = 7t, xeD, te]0,T], 



(2.3) 



w{t,x) =0, X e 9D, t e]0,T], 

Ji3w(t,x)dx = f(t), te[0,T], 
^w(0,x) = wo(x), xeD. 




We follow essentially the same approach of [ ], with additional care on the analysis of the 
initial datum and on the dependence of the solution on the various parameters of the problem. 
In the sequel we will employ a spectral (spatial) approximation using the L^(D)-orthonormal 
eigenfimctions (ek)keN of the Laplace operator, 

xeD, 
X e 9D. 

Define 

(3j:=(l,ei), jeN, 

where ( • , • ] denotes the (D ) scalar product and 1 is the function defined on D such that 1 (x) — 
1 a. e.. Without loss of generality from now on we assume that |D|, the Lebesgue measure of D, 
is equal to one. Clearly, 

oo oo 

l = ^(3jej(x) and ^ (3^ = [D| = 1. 

k=l k=l 

A special role is played by the pipe's flux carrier, i. e., by the function O (which belongs to 
H^(D) n HJ(D) under the smoothness assumptions on D) defined as solution of the following 
Poisson problem 




xeD, 
X e 9D. 



We define the following quantities 



O dx = 

D 



iVOrdx > and ri;^ 

D 



lordx > 0, 

D 
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which clearly depend only on D. Observe that the function O is enough in the stationary case 
to construct Poiseuille-tjrpe flows, since in that case the problem for the flux and that for the 
pressure completely decouple. On the other hand, in the time-dependent case the situation if 
more complex, since both unknown depend also on the time. 

To work with in our problem with a general viscosity we need the scaled version of the flux 
carrier cp :— v^^O which solves 



(2.4) 

and such that 




cp dx = V 



|V(ppdx= ^ >0 
D y 



xeD, 
X e 9D, 

and 



i(prdx = 



-v/2' 



xs 



We start our analysis by proving the following result. 

Proposition 2.7. Given f e H^(0,T), assume that^ wo(x) = 
solution (w, 7t] of (2.3) such that 

we C(0,T;Hj(D))nHi(0,T;L2(D))nL2(0,T;H2(D)), 

7teL2(0,T), 

satisfying the following estimate 



^'P|£lf(o). Then, there exists a unique 



v||Vw(t)|| 



(2.5) 



Aw(s) 



+ 



ds 



7T S 



ds 



^ C 



((l+v2)|f(s)P + (l+v)|f'(s)P) ds, 



with a constant C depending only on D (and in particular independent ofJ). 

Proof. We start by constructing, with the Faedo-Galerkin method, a global unique approximate 
solution in V-m — Span(ei, . . . , era). The first step is to approximate the initial condition. Let 
(p be the function introduced in (2.4) and write <p = ^^Li VkGic/ where the series converges in 
H^(D) and cpic — (cp, eic). Hence, the projection of cp over Vm, is given by 

m 

k=l 

In order to satisfy the flux condition also at time t = we set 

w^{0,x) -p— r--T-- 

Observe that, for large enough ra G N, the approximate initial datum is well-defined. In fact, 
Pm^P — > (p in L^(D) and since |D| < +oo then Ptr9 — > (p in L^(D). Since Pm^ ^ So ^ ~ 
'^^^Xo > 0' there exists mo G N such that Pm^ 7^ for all m ^ rao. Moreover, w™-[0, x) — > 
Wo(x) in H^(D), as m — > +oo. 



We write Galerkin approximate functions 

w'^(t,xl = 



k=l 



^The initial condition here is chosen in such a way that the compatibility conditions on the flux at time t = are 
satisfied. 
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and we look for a couple (w^, tt^) such that 



(2.6) 



— (w-,ej)+v(Vw"^,Vej) =7t™(l,ej) for j = 1, . . . , m, 



and : (0,T) (0,T) chosen so that the flux condition 



(2.7) 



w"^(t,x) dx = f{t) Vte(0,T), 



is satisfied. The equality is meaningful since f is a. e. equal to a continuous function. In terms 
of Galerkin coefficients (Cj^Ji-jj^m we have for the initial condition that 

while the system of ordinary differential equations reads as 

Ac-(t) + vAjC-(t) = 7T-(t)Pj for j = 1, . . . ,m, 
and the solution can be written as follows 



cr(t) = cr(o)e- 



-vAit 



7t"^(s)e 



-vAi(t-s) 



ds for j = 1, . 



, m. 



To find the equation satisfied by tt"^ we multiply the latter equality by (3j and sum over j 
1, . . . , m to get 



j=0 j=0 



7t"Ms e 



-vAj (t-s) 



ds. 



Finally, to obtain an integral equation for tt"^, we differentiate with respect to time deducing 



j=0 j=0 

and this yields the following Volterm integral equation of the second type 



ds , 



(2.8) 



t m ;x.(d2 



^-vAit 



l(3F 



where |(3p :— ^J^^ (3^. For any fixed m G N the kernel of the integral equation (2.8) is bounded 
for all ^ s ^ t. This is enough to infer that if f G H^(0, T), then there exists a unique tt™ G 
(0, T) satisfying the integral equation (2.8) and such that 

||7t"^||l2(o.t) C^(v)||f|lHMo,T), 

for a constant Ctti("v) possibly depending on m and also on v. Especially the dependence on 
m is crucial, since we consider the problem at fixed viscosity, while we need uniform estimates 
in m G N to employ the Galerkin method. In particular, the uniform estimate does not follow 

directly since the series defining the kernel for s = t, that is(^J^i|3^) ^ Y.'v=i^^^v does not 
converge for m +oo, see [ ] for further details. 

We need to find the a priori estimate in a different way, but observe that, once we have con- 
structed TT™^, we can use it as a given external force in the equation for the velocity (2.6). By using 
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w™- as test function, we obtain (with the Schwarz inequality and by using (2.7)) that 



Vw 



m||2 



7t""(s)(l,w"^{s))ds 



7t"^(s)f(s) ds 





1 



^ 



|7T"^(slMds. 



Observe also that wj^ wq in L^(D) and that ||wj^|p 2||wo|p for large enough ra and in 
addition 

||Vw-||2 <; 2||Vwo|p _||f||^,,„^^,||Vcp|p 2||f||2„,,o.T-,. 

AO 

This gives the first a-priori estimate showing that, for all m ^ mo there exists a unique solution 
w"^ e L°°(0, T; L2)nL2(0, T; H^). Since the bounds we obtain on the solution are not independent 
of m, they cannot be used to make directly the Galerkrn method to work. We can prove even 
more regularity on w"^ by standard estimates. In fact, by using as test function in the system 
satisfied by w"^ we get 

llwrWiP + ^^l|Vw-(t)||2 = 7r-(t]{l,wr{t]) = 7r-(t)f'(t). 



Hence, an application of the Schwarz inequality gives 



(2.9) _||Vw-(t)||2 + 



IwlT'lpds «:a/||Vwo|P 



2e 



If'fslPds. 



By using — Aw"^ as test function in the equation satisfied by w"^ we also obtain 

'||Aw"^(s)|pds^ ||Vwo|p + -^ 

> ^'^ 



(2.10) 



i||Vw-(t)|p + | 



ds. 



These estimates are enough to construct, for each fixed m e N, a unique solution (w"^, tt"^), 
which is smooth, say 14^(0, T; 1^(0)) n 1^(0, T; 14^(0)). Again the presence of tt"^ in the right- 
hand side prevents from imiformity in ra. 

This issue is solved by using a special test function and a couple of nested a-priori estimates 
in the following lemma. 

Lemma 2.8. There exists Mq G N (larger or equal than roo) such that for all ra > Mq it holds for all 
te [0,1] 



v||Vw"^(t)||2+v2 



(2.11) 



I Aw" 



'ds 



7T" 



ds 



s$ C 



((l + Y2)|f{s)P + (l+v)|f'(s]P)ds, 



with a constant C depending only on D (hence independent ofm and of J). 

Proof. We observe that the fimction Pm V, where cp is defined in (2.4), turns out to be a legitimate 
test function for the Galerkin system. With integration by parts, we obtain 

By adding to both sides of the previous equality the quantity — v(w"^, Acp) — (w^, 1) — f (t] we 
obtain 

(W^, Pm(p) + f = 7t"^(l, cp) + 7t^(l, PmCp - cp) + ^'(W"^, A(P^(p - (p)). 
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The last term from the right-hand side vanishes since e Vm, while A[?rnV — cp) € V^. 
Hence squaring the latter equality (remember that [l,(p) — v^^Xo) obtain 



|7T"^{t)P ^ 3(||w-(t] 



|f(t)P + |7r-{t)p||P^cp-cpf). 



Next, since Pm is a projection operator it follows that there exists Mq G N such that UPmip — 
"pIP < Xo/(6v^) for all m ^ Mq. Consequently, we can absorb in the left-hand side the term 
involving tt™ from the right-hand side. Consequently, after integration over (0, T), we get 



(2.12) 



7T" 



(sir ds c 



(II 



-v2|f(s)P) ds Vm^Mo, 



with a constant C depending only on D {via xo and rjo). Consider again (2.9) and use now (2.12) 
with e — 1/C. It is possible to absorb the term involving n from the right-hand side, obtaining 

1 rt 



< vIlVwol 



C 



'ds 



v2|f(s)P + |f'(s)Pds 



^ C 



{l + v2)|f(s)P + (l+v)|f'{s)Pds 



This shows a uniform bound on||w]^||L2(o,T)- Hence, coming back again to (2 . 1 2) we also obtain 
that 



7t 



(sir ds c 



(1+Y2)|f(sl|2 + (1 + V)|f'(s)|2ds 



Next, by using the bound obtained in (2.10) we also get 



|Aw""(sll|2 < C 



(l+v2l|f(s)|2 + (l+A/)|f'{s)Pds 



□ 



By collecting all the estimates and with Young's inequality we obtain (2.11). 

With the above lemma we can conclude the proof of the existence result, since uniform 
bounds imply that there exists a couple (w, tt) and a sub-sequence {rakjkeN such that 

inHi(0,T;L2(D}}nL2(0,T;H2(D)), 



w 



7T 



^ W 
^ W 
7T 



inL°°(0,T;Hj(D)), 
inL2(0,T). 



The problem is linear and this implies that (w, tt) is a distributional solution of (2.3) with the 
requested regularity. Uniqueness follows again from linearity of the problem. □ 

Since the estimates on the norm of the solution are independent of T (they depend just on 
D and on v) the same argument shows that if f € H^(0, +oo), one can study the problem with 
arbitrary T > and then a unique solution (with the same regularity) exists in [0, +oo). More 
generally we have also the following result on the whole real line, which is obtained by letting 
the initial time to go to —oo. 

Corollary 2.9. Let he given f e (R). Then, there exists a unique solution (w, Tt) of (2.3) defined for 
flZZ t e R such that 

w e Co(R; Hj(D)] n H\R- L^ID)] n \3[R-H^[m, 
Tte l2(R), 



ON LERAY'S PROBLEM FOR ALMOST PERIODIC FLOWS 



11 



with the same bounds as before (here Co denotes the subspace of continuous functions vanishing at 
infinity). 

For our purposes of studying almost periodic solutions it is important to show that one has 
a global solution, with uniformly bounded gradients, also if the force f is not in (0, +00), but 
just in H^ig^(0, +00), where 



|f||Hljo,+oo) := sup( 



■t+1 



|f(s)P + |f'(s)pds 



In particular, the following result will be crucial for the rest of the paper. 

Proposition 2.10. Let be given f e Hj,^^(0, +00), then the unique solution of (2.3) exists for all positive 
time and it satisfies 

w e Cb(0, +00; Hi(D)) n HliJO, +00; L2(D)) n ^,^^0, +00; H^{D)), 
with the estimate (2.19) in terms of the data. 

Proof. Generally this result is straightforward in presence of a standard parabolic problem. 
Since here we deal essentially with an inverse problem, we give a detailed proof, which is 
nevertheless obtained adapting the usual techniques typical of almost periodic solutions, see 
e. g. Amerio and Prouse [1]. Observe that the estimate (2.11) does not give a direct control of 
supo<t<T II Vw(t)|| since the bound depends on ||f||Hi(o,T) arid consequently the H^-norm of w 
may become unbounded when T +oo, if f ^ H^(R). We first prove that ||Vw|| e L°°(0, +oo). 
First observe that Proposition 2.7 imply that there exists a unique solution 

w e C(0, +oo; Hi(D)) n Hl,[0, +oo; LHD]] D L^^JO, +oo; H^D)), 
7tG Lf„JO,+oo), 

hence the following calculations will be completely justified. By using the function cp from (2.4) 
as test function we obtain with integration by parts the following identity 

(wt, cp) — v(Aw, cp) — (wt, cp) — (w, vAcp) = 7t(1, cp]. 

Hence, by recalling the definition of cp and the flux condition we get 

^7t(t) = {Wt(t,x),(pW)+f(t) 
■V 

and consequently taking the square and integrating over the [£,, t] (for all couples £, ^ x) we 
obtain that there exists C depending only on D such that 



(2.13) 



|7t(s)P ds^C 



wt(s)|r +v1f(s)r ds. 



Next, we test the equation satisfied by (w, tt) by w and v ^Wt, and by recalling that Wt dx — 
( JjjWdx)' — f'(t) we obtain (for any e,ri > 0) the following differential inequalities, a.e. 

t e (0,+oo) 

i|Hp + ^||Vwf ^||7tp + J-|fp, 
i||Vw|P + i||w.|r^AH. + _l.,f^,.. 
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By choosing e — [2yC] ^ and r| = (2C) ^, where C is the constant appearing in (2.13) and by 
integrating over an arbitrary interval [£,, t] we get 



W T 



HI 



(2.14) 



(v||Vw(s)||2 + -||wt(s)||2)ds 



^'(2 + C)|f(s)|2 



-|f'(s)Pds. 



Hence, by using the Poincare inequality and dropping the non-negative term v 
finally obtain that there exist Ci, C2 > and depending only on D such that 



wt we 



(2.15) 



IIwItJIIhi +CiV 



\\Ms]\\hi ds < i||w{t)||^i +C2 



v|f(s)|2 



1 



|f'(s)P ds. 



Suppose now that for a given t e [0, +00 [ it holds that 
(2.16) ||w(t)|p„, ^ |lw(t + l)fH, 

and rewrite (2.15) in the interval [t, t + 1] as follows: 



(2.17) ||w(t + l)||2,i - \\yv{t)\\f^, + CiV 



w s 



HI 



ds ^ C2 



t+i 



v|f(s)P 



1 



|f'(s)Pds. 



Since by hypothesis (2.16) the first term is non-negative we obtain in particular that 
(2.18) CiV 



t+i / ^ 

||w(s)||^i dss;c2 + - 



1^^ 



Next, by using the estimate (2.18) and the same argument we get, for each couple t ^ £, with 
T, £, e [t,t + l],that 



■T 1 

||w^(s)||^i ds + C2(v + - 
^2c2(v+i)||f||2 . 

"V uloc 



Since ||w||^i is a continuous function we can fix f, e \t,t + 1] such that 



Hi 



mm ||w[sJ||hi. 
t<s<t+l 



We have then, using (2.16) and the definition of 

\w[t+i]\\f,,~\\MmH^ + \\Mm 



|lw(t)f„, ^|lw(t + l)||^, 



2 

HI 



IIM^)IIhi 



■t+1 



|w(^)||^i ds 



|w(t+l)| 



2 

HI 



+ 



w s 



I HI 



ds. 



Finally, from (2.17) and (2.18) we obtain 



Hence, for each t e (0, +00) such that (2.16) holds true we have that ||w(t)||^i is bounded 
uniformly. On the contrary, if 



iw{t)rH, >iiw(t+i)ii 



Hi' 
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we can repeat the same argument on the "window" [t — 1 , t] . In this way it is clear that if we are 
not able to find an interval [t — n— l,t — n] with n < t — 1 such that 

||w(t-n-l)||^i ^ ||w(t-n)||^i 

we obtain, for some n e N such that ^ t — n < 1, the inequalities 



||w(t-n) 



HI 



> 



> Mt] 



hence that 



||w(t)||Hi ^ sup ||w(t)||^ii 
O^t^l 



The latter is bounded simply by the local existence result. This finally shows that, since || wq || hi 
depends itself on the norm of f in H^[|j^(R), for all t e R+, 

rt+l 



w(t) 



||w(s)||f^l ds C v + 1 



1 



Next, by going back to (2.14) we obtain an estimate on ||wt||L2(-|- t+i)' which implies by (2.13) 
the corresponding estimate for |7t|^2(t t+i)- Fu^ally/ by comparison we get also an estimate for 
II Aw||L2(t and by collecting all inequalities we finally get 



(2.19) 



sup 

t>0 



v||Vw(t)|| 



t+i 



(v2||Aw{s)|p + ||wt(s) 



7T S 



ds 



1 



1. 



ending the proof of the Proposition. □ 
We can finally construct a solution over the whole real line and we have the following result. 
Proposition 2.11. Let be given f e Hj^,g^(R), with 

1/2 



|f|lH,V,(R) - sup 

t fc K, 



t+i 



|f(s)|2 + |f'(s)pds 



then the unique solution of (2.2) exists for all times and it satisfies 

w e Cb(R; Hi) n Hi,„,(R; L^) n L2„^(R; H^), 

Tte l2,„^(R), 

satisfying the estimate (2.1). 

Proof. We start by observing that since f G HJjq^(R), then it follows that f e Cb(R) (again by 
a.e. identification). It follows directly that f e C(R), but the control of the maximum of |f| is 
obtained as follows. We claim that sup^^p^ |f(x)| ^ 2||f||Hij ■ In fact, for each couple of points 
X, y € R such that |x — y | ^ 1 it follows that 

py 

f'fs] ds 



|fW-f(y)| 



^ llfll 



Suppose now per absurdum that there exists xq G R such that |f(xo]| > 2||f||H^ • The previous 
inequality implies that 

|f(x)| > mnl,^^ forallxe [xo,xo + l], 

hence the contradiction 



xo + l 



lf(s)pds>||ff^, 
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This proves that the bound on |f (x)| is true for all x G R. 

The proof of Proposition 2.11 is then obtained by using the previous results from Proposi- 
tion 2.10 to solve the following family of problems parametrized by n G N, 



Wn(t,x) = 0, 

Jo Wn(t,x)dx = f(t) 
[w^(-n,x) = ^f(-n), 



X G D, t > -n, 
X G 9D, t > — n, 
t > — n, 
X G D, t = — n. 



The same arguments as before imply that 

w^GL2i„J-n,+oo;H2(D))nHii„J-n,+oo;L2(D)), 
Vwn G L°°(-n, +oo;L2(D]], 

with bounds independent of n G N (the dependence on v is the same as in (2.5)). By defining 
the extended functions Wn(t, x) and 7Tn(t, x) on the whole real line as 



Wn(t,x) <^ vy(^) 



Xa 



f(-n) 



t ^ -n, 
t -n, 



and 7Tn(t) := 



7tn(t) 

vf(-n) 



t ^ -n, 
t < -n, 



we can extract a sub-sequence (relabeled) as (wn (t, x] , 7tn (t, x) ) such that 

- w in L^i^^ (R; H2 ( D ) ) n H^i^, (R; ( D ) ) , 



VWn ^ Vw 



inL°°(R; L^ID)), 



inl-^i^JR) 



It is easy to see that (w, tt) is a distributional solution to (2.3). Since the problem is linear this is 
the unique solution. □ 

We finally prove the result in Stepanov space of almost periodic functions. 

Proof of Theorem 2.2. Since the function f is almost periodic, for each sequence [rm) C R we can 
find a sub-sequence {Ttr^} C R and a function f such that 



sup 
teR J 



|f(T + r, 



If'fT + r, 



0. 



By using a standard argument by contradiction (see e. g. Foias and Zaidman and Foias and 
Prodi [ ]) we assume per absurdum that w and n are not almost periodic, hence that there exist a 
sequence {Km.} C R, a function f : R ^ R such that 



sup 
teR J 



(If (T + h„,) - f (t)P + If '(T + h^) - f '{t)|2) dx 



0, 



a constant Sq > 0, and three sequences {tp}, {hrap}, {Hn^J such that for all p G N, 



sup 
teR 



(v^||Aw(T-hh^ 



Aw(t + hn 



+ |7TfT-hh^ 



7t(T + h„ )P) dT^ 6o. 
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In addition, there are (by eventually relabeling the sequences) two real functions f i and f2 such 
that 



sup 

tGR J 



sup 

t6R 



t 

■t+1 



(If (T + tp + J - fiWP + If '(T + tp + J - (t)P) dT 0, 
(If (T + tp + J - f2(T)P + |f'(T + tp + J - f^(T)l2) dX 0. 



It is clear that f = f i — ^2- We consider now the problem (2.3) with the two fluxes 

Fip(t) :=f{t + tp +h^J, 

F2p(t) :=f{t + tp+hnj, 

and the corresponding solutions 

wip (t, x) w(t + tp + h^p , x) , Ttip (t, x) n{t + tp + h^p ) , 

W2p(t,x) :=w(t + tp +h,np,x), 7t2p(t,x) := 7t{t + tp +hnp). 

Observe that Fip(t) — > f(t) and F2p(t) f(t) in H^[^^(R) as p oo. Hence, by passing to the 
limit as p — > +00 we construct two solutions (wi, Tti) and (w2, 7T2) corresponding to the same 
flux f . In particular, 

Wip — >• Wi and W2p W2, 

in Cb(R.; Hj(D)), but also in the topologies given in the statement of Proposition 2.11. Hence, 
we have in particular 



60 < 



v2||Aw(t + h^p^ ] - Aw(t + Knp Jf + |7r(T + h^p^ ) - 7t(t+ h^p^ )P dx : 



■■PS - 

fl 





V^ll Awip, (X) - Aw2p, {X)|p + |7Tip^ (X) - 7t2p, (x)P dX 



because the problem is linear and the two solutions wi and W2 corresponding to the same flux, 
coincide. This proves the almost periodicity of w and the same argument applied to Wt ends 
the proof of the result. □ 

2.3. The full nonlinear problem in S^{R). We now finally consider the Navier-Stokes equa- 
tions and we look for solutions of the following problem 

'atU-vAu+{u- V)u + Vp =0, xeO, teR, 

V - u = 0, X e O, t e R, 

u = 0, xe90, teR, 



with 



sup||u{t) -Wi(t)||Hi(Oi) Co, 
teR 



for i = 1, 2. Observe that this constraint implies (see [ , § 7] and [ , VI.I]) that 

lim ||u(t) -Wi(t]||Hi/2(Di)) = 0, 

uniformly in t. 

As a preliminary remark we recall that it is well-known (see for example [ , VI], Amick [ ]) 
that due to the particular shape of the unbounded domain O the Poincare inequality holds true 

(2.20) 3cp>0 ||u||l2(o) <Cp||Vu||l2(o) VueHjlO) 
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and, if the boundary is smooth enough (say of class C^'^), then the following estimate for the 
Stokes operator holds true 

|1u|1h2(o) cl|PAuj|L2(o), u e h2(0) n HliO). 

We now take advantage of the results on the linear case of the previous section, but first we 
need to "glue" the basic flows constructed in the two pipes. In this part of the paper we do not 
claim any originality and we use a classical approach. We denote by z € R"*" the axial coordinate 
in both cylinders and we define the "truncated pipes" 

0[:={(x,z) e Oil z<r} i=l,2, 

and we also define the truncated domain 

:= OoUOIUO^. 

We first define a field (which is as smooth as wt) defined on O and such that is equal to Wt in 
the sets Oi\Oj. This extension is obtained by freezing the time variable and by gluing together 
the functions wt by cut-off functions ^i{z) e C°°(R^) depending just on the axial coordinate z, 
and such that 

fi xeOAO{, 
= \o X e oY\ 

Next observe that, 

2 

(2.21) Vo := ]^<t>iWie Hii„,(R; L\0)) n L^^^JR; HHO)) 



i=l 

since the mapping (wi,W2) Vq is bi-linear and the extension does not involve the time- 
variable. 

Remark 2.12. We also observe that if we are in a different functional framework the same pro- 
cedure can be applied because the properties of the extension with respect to the time variable 
are the same of those of the basic flows (w^, 7t^). 

This is not exactly the required extension, since the function Vq is not divergence-free. To this 
end one has to use the Bogovskii formula to solve the linear problem in the boimded domain 
O^: Find S Vo such that 

fv- (SVo) =-V-Vo, xeO^ 
|3Vo = 0, xeao^. 

Since V • Vq G Hg(O^) and the compatibility condition Jqi ^ • = is satisfied the problem 
has a solution !B Vq e (O^ ). Denoting again by IB Vq the null extension of O^, we finally set 

(2.22) w:=Vo + ^Vq. 

The regularity of the solution of the divergence equation and the previous argument shows that 
if Wt € Hii<,JR;L2(0t]] nL^„,(R;H2(0|)),then 

w=Wt in Ot\Oj 

and 

2 

||w||Hi(R;L2(0))nL2(R;H2(0)) < C ^ ||Wi || hi^JR;L2 (Ot ))nLS„JR;H2 (Ot )) > 

i=l 

for some c = c(0o, Oi, O2). We have finally the following result. 
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Lemma 2.13. Let (wi,,7Ti) be solutions of the basic flow with the regularity of Theorem 2.2. Then the 
extended function w is Stepanov almost periodic and 

l|w||L'»(R;Hi(0)nH,V(R;L2(0))nL2„^(R;H2(0)) 
2 

^ ||Wi||Loo(R.Hi(o)nH„V-(R;L2(Oi))nL2„JR;H2(Oi))- 

1=1 

With this result, we look now for solutions in the form 

U = U + W, p=P + 7T. 

We have to find U e L^^JR; l-^[0)) D L^j^JR; Hi(0)) solving 



(2.23) 



9tU - vAU + (U • V) U + (U • V) w 

vu = o, 

U = 0, 
[U(0,x) = -w(0,x). 



(w- V)U + VP = F, 



on 90, 



where 



(2.24) 



Fft,x,z] 



-(9tw(t,x) - vAw(t,x) + (w(t,x) • V)w(t,x)) + ^ — (z4)i(z]7ti(t)] 



i=l 



From the results of the previous section we can infer that F e S^(R; L^(0)), hence we can now 
use standard results to show existence of an almost periodic solution U. The support of F is 
contained in the bounded subset 0^, hence we can use the standard variational techniques to 
show existence of a unique solution, provided that the viscosity is large enough. In particular 
the only property that we need to check is that (w(t,x) • V)w(t, x) is almost periodic. This 
follows since if we take the V{0] -norm of F, this is equal to the L^(Oi)-norm. All other terms 
in the summation are clearly relatively compact in L^jj^^(R], the nonlinear one can be estimated 
as 

rt+1 f 2 pt+l 

w(t,x) • V)w(t,x) ^ sup|lVw(t)||2,(Q^) ||Aw(s)||^2(o,)ds. 

Jt JOi t Jt 

This proves that F e S^(R; L^(0)) and, by using this expression and from the estimate (2.1) 
(w^, 7t^), we obtain that 



on 



IIfII'l^, <c 



y 



If|l^.^ 



In particular, as in the previous section, we need to show that there exists a solution U € 
L^[g^(0, +oo; H^(0)), and this will follow by using that the right-hand side is in l-^^lg^[0, +oo; L^(0)). 
The existence of a local solution Lj^ocf*^' +°°; Hj{0)) follows by standard arguments. In partic- 
ular, one has to perform a truncation in the space variables (which is possible since F is of 
compact support, see again [ ]) and the usual a-priori estimates. In addition to these rather 
standard results, we now prove that the solution is strong (provided that the viscosity is large 
enough, which is nevertheless needed also for weak solutions) and that the solution is uni- 
formly bounded in H^(0) and L^2j,j,(R; H^(0)). We present just the a priori estimates, which can 
be justified by the usual Galerkin method and truncation of the domain. 
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Proof of Theorem 2.4. We multiply (2.23) by — PAU and we carefully treat the various terms. First 
note that 

fk+l 



Oi 



(PAU] • (U- V)U = ^ 



(U- V)UPAUdxdz 



Di 



< C||U||l4((^ j^+l)xDi)ll VU||l4((^ k+ljxDi) l|AU||L2((k,k+l)xDi) 



^ C||U||hi(o,)I|VU||hi(o,)II^LI||l2(o,), 

since the constant of the Sobolev embedding c L"' are uniformly boimded in each of the 
strips (k, k + 1) x Dt, see [ , Lemma 2.1]. Moreover the standard regularity theory for the 
Stokes operator shows also that 

(PAU) • (U • V)U ^ C||U||hmo,)I|PAU||2,(o^,. 

The term Jq (U • V) w is then estimated as follows, by using the same splitting into slices of 
width 1 in the z-direction and the Sobolev embedding c L°°, 

(U.V)wPAUs^C||w||l2(o,)I|PAU||2.,o^), 

The other term is estimated in the same way, and the contribution from the domain Oq is han- 
dled with standard tools. Hence adding together the three terms, we finally we arrive at the 
differential inequality 

A||vU|p + (v - Ci(||VU|| + ||Vw||))||PAU||2 < 

We fix the viscosity large enough, so that 

V- Ci(||VU(0)|| -||Vw(0)||) =A'-2Ci||Vw(0)|| > 0. 

and this is possible since ||Vw(0)|| = ^^|f(0)| sC vl|f|lH;i„- Since ||F(t)|| e L^j^JR) this is 
enough to show existence of a local unique solution, in a time interval [0, 5 [ for some positive 5. 

The next step is to show that, under the same assumptions, the solution is global. This is 
rather standard, since if we have the uniform bound on || Vw|p coming from Proposition 2.11, 



sup||Vw(t)f < c(v + l + — )||f||^ 



tGR 



V^ 



JR.) 



we can fix — Vq (f , O) large enough such that 



ifv>Vo then v - Ci||Vw(t)|| ^ 



t e R. 



For such y we are reduced to solve (possibly redefining the constants) 



^livuir 



ci( 2 -||VU|| 



IVUII^ < ^||F||2 



Hence, if we define Z(t) as the solution of the following Cauchy problem 

.2 25) ^Z' + Cl(f-^/Z)Z=f||F|P 

\Z(0) = ||VU(0)||2 

we have that ||VU(t)|p ^ Z(t). We employ now a fixed point argument in the space of function 
which are continuous and boimded over [0, +00] to show that Z is well defined and bounded in 
the same interval. Let be given z G Cb (0, +00) such that 
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Solve now the problem 

fz' + Ci(^-^^)z=^|lFf. 
\z(0] =Z(0) 

Since ^ — CiVz. ^ ^ and both z(0) and -^HFH^ are non-negative, it follows that z(t) ^ 0. Hence 
z satisfies 

(2.26) z' + ^z^^||F|P. 



The same argument employed in the proof of Proposition 2.10 shows that 
Then, if v is large enough such that 

c(v + l + 4)l|f|lHi 

Consider now the map O : z i-^ z and given t > define 

Kt::= jze C°([0,t]) : z ^ - 



Clearly we have O(Kt^) C K^^ and the map is relatively compact by Ascoli-Arzela theorem since 
z is Lipschitz continuous as solution of (2.26). Hence O as a unique fixed point which is a 
solution to (2.25). Since x > is arbitrary this proves that Z exists on the whole interval [0, +oo). 
A standard comparison argument shows that || VU(t)||^ ^ Z(t) ^ ^ for all positive t. 

This estimate implies, by using standard argument well established for the Navier-Stokes 
equations, that there exists a solution 

U e Cb(0, +oo; Hj(0)) n L^i^jO, +oo; W^O]] n HlJO, +oo; 1^(0)), 

and due to the regularity proved the solution is unique. 

Next, one can construct a global solution, by solving the following family of problems in 
[— n, oo) X O, 

f 9tUn - vAUn + (Un • V]Un + (Un • V)W + (W • V) Un + VPrx = F, 

V • = 0, 
Un = 0, on 90, 
^Un(— ti, x) = — w(— n, x) in O. 



Again by prolongation we define a velocity on the whole real line by 

Un(t,x) = 



Un{t,x) t^-n, 
Uf— n, xl t < — n. 



We can show that 11^ (t, x) converges to a solution U such that 

U e Cb(R; Hj(0)) n L^^^JR; H^fO)) n H^JR; \}[0)). 

To end the proof we need to show that if the external force is almost periodic, then U is almost 
periodic too. For this result we need a result of "asymptotic equivalence," which is obtained as 
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follows. Let us suppose the we have two solutions of (2.23) on the interval [to, +00) correspond- 
ing to different initial data, but to the same external force. Then the difference U := Ui — U2 
satisfies 

( 9tU - vAU + (Ui • V) Ui - (U2 • V) U2 + (U • V) w + (w • V) U + VP = 0, 

v-u = o, 

U = 0, on 90, 

By multiplication by U and usual integration by parts, we get 
1 d 



2 dt 



2 + v||VU||2 <; 



(U-V)(w + U2)Udx. 



By using the same techniques employed before to handle unbounded domains we show that 

i A||U||2 + C(||VW|| + ||VU2||))||VU||2 sc: 0. 

The uniform bounds on w and U2 in H J ( O ) imply that for large enough viscosity one has 



A||uf +vCp||U|p<CO. 



This finally shows that 



|U|l^^||U(to)re 



2 -vCp(t-to) 



t ^ to. 



The same argument employed in the proof of Theorem 2.2 can be employed and, for each se- 
quence {Ttti} C R we can find a sub-sequence {tra^ } C R and a function F such that 

rt+l 



sup 

t6R J 



|F(T + r^J-FW||2dT'^^°°0. 



Assume by contradiction that U is not almost periodic, hence that there exist a sequence {hm} C 
R and a function f : R R such that 

rt+l 



sup 
teR J 



|F(T + h^)-FW|PdT"^^°° 



0, 



and a constant 60 > and three sequences {tp}, {hrap} and {h^p} such that for all p G N, 

■tp+i 

sup ||U{T + h„, J-U(T + hnp)|pdT^5o. 



There exists two real functions Fi and F2 such that 

rt+l 



sup 

tGR J 



F(T + t„+h„ l-FiWII^dT'"^ 



sup 

tGR J 



t 

■t+1 



|F(T + tp 



F2(T)|pdT'^^°°0. 



It follows in a standard way that F = Fi = F2 and we consider now the problem (2.23) with the 
two forces 

Fip(t,x) := F(t + tp +hn^p,x] and F2p (t, x) := F(t + tp + h^^ , x], 
and the corresponding solutions 

Uip{t,x] :=U(t + tp +hn^p,x) and U2p (t, x) U(t + tp + , x). 
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By passing to the limit as p — >^ +00 we construct two solutions Ui and U2 corresponding to the 
same force F. In particular 

wip ^- wi in Cioc(R; L^), W2p ^- W2 in Cioc(R; L^)- 

In addition we have 



60 < 



||U(T + hn^^J-U(T + hn,J||2dT: 



rl 



||UipjT)-U2p3WfdT 



|Ui(T)-U2(T)||2dT. 



On the other hand the asymptotic equivalence implies that Ui = U2, since for all t e R, 

||Ui(t) -U2(t)f < ||Ui(to) -U2(to)||2e-^c-(t-t°) t > to, 

and letting to —00 we obtain a contradiction. 

To conclude the proof we show that the function U € S^(R; H''(0)) for all ^ s < 2. In 
fact, take a sequence {r^} we can find a sub-sequence {rm.^} c R and a function U such that the 
sequence {U(t + )} satisfies 

rt+l 



sup 
teR J 



||U(T + r^J-U(T)f dT'^^°°0. 



Hence the sequence {lI(T + rTn,k)}is a Cauchy sequence in L^q(.(R; L^(0)), that is for every e > 
there is N e N such that 

-t+i 

sup ||U(T + r^J-U(T + r^J||2dT< e p,s>N. 
teR Jt 

Since U e H^^JR; ^{0)) n L^oc^^*'; by classical interpolation it follows that 

U* e C(R;m(0,l;H2-2'^(O))), < cr < 1. 

Hence we obtain that 

rt+l 

||U(T + r,^ J - U(T + Jll?^. dT ^ 

||U(T + r^ J - U(T + r^3 ) ||^7fo) l|U(T + J - U(t + r^ J || (q) dx, 

and by using Holder inequality, 
rt+l 



sup 
teR J 



< sup 
teR J 



||U(t + J - U(t + J||2,3 dx < 



t 

t+i 



||U(T + r^J-U(T + r^J||(.(0) 



||U(t + J - U(t + t^Mhho) dT 



Since 



sup 
teR J 



l|U(T + r^J||^.(0) + ||U(T + r^J||^.(o) dx < C, 



it follows that the sequence {U(T-l-rTn,k )} is a Cauchy sequence in L^jj^^(R; H^(0)) as well, ending 
the proof. □ 



Remark 2.14. The result with f e (R) follows in the same, even simpler, way. 
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3. Leray's problem in the framework of Besicovitch a. p. functions 

In this section we discuss the same problem in a more general setting, and first we recall 
some definitions on almost periodic solutions. 

3.1. Generalities on almost periodic functions. In the literature there are different definitions 
of almost periodic functions and we need now to explain the precise setting we are using. We refer 
mainly to [ , ch. I] for further details and references. Let Trig(R) be the set of all trigonometric 
polynomials, that is, u e Trig(R) if there exist n £ N, £,i , . . . , f,n £ R and Ui , . . . , £ C such 
that 

n 

u(x) = ^Uke'^"", xeR. 

k=l 

Next, a set A C R is relatively dense if there exists L > such that each interval of length L 
contains an element of the set A. 

Definition 3.1 (Bohr). A uniformly almost periodic function (UAP(R)) is a continuous function 
f : R R such that there is a relatively-dense set of e-almost-periods. That is for all e > 0, there 
exist translations > of the variable t such that 

|f(t + T,)-f(t)| SC £. 

It is easy to see that all trigonometric polynomials are almost periodic according to the previ- 
ous definition. Let UAP(R) be the set of all uniformly almost periodic functions. Then UAP(R) 
coincides with the closure of Trig(R) with respect to the sup-norm \\ ■ ||l°=. Alternatively, as re- 
called in Section 2, a function f G U AP (R) if the set {f (t + • ) : t G R} of translates of f is relatively 
compact in C(R). A more general notion of almost periodicity was introduced by Stepanov in 
1925. To this end for p ^ I and r > 0, define the norm 



f||sp,T := sup( - 
teR 



Then, the space SP(R) is the closure of Trig(R) with respect to the norm || • ||sp.r above. Notice 
also that while the norm depends on r, the topology is independent of the value of r, hence we 
re-obtain the definition used in Section 2 (Cf. [ ]). 

The definition was later extended by Weyl in 1927 by considering the closure WT'{R) of 
trigonometric pol5momials with respect to the semi-norm 

||f||wp lim |if|!sp.r- 

r— >oo 

Finally, Besicovitch [5] defined the space 23^ (R) as the closure of Trig(R) with respect to the 
semi-norm 



( 1 



|f(s)P ds 

R 



||f||sp :=liinsup( — 

Notice that one can have II f II sp = even though f ^ 0. For example this happens if f is in L'' (R) 
with q > p or f G L°°(R) and |f (x)| — >^ as |x| oo. One has the following strict inclusions 

UAP(R) c SP(R) C WT(R) c ST'(R), for any p G]l,-hoo[, 

(with obvious inclusions with different values of p). It turns out that the spaces (R] gf Besi- 
covitch almost periodic functions are among the "largest possible" compatible with the treat- 
ment of partial differential equations as we shall see in Proposition 3.2. Let us focus on the case 
p = 2, since 23 ■^(R) has an Hilbert structure. Given f G Lj^^^^fR), define 



M(f) := lim sup 

R-s-oo 



1 '■'^ 



f (t) dt. 

-R 
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The mean operator M(f ) is defined as the above quantity when the limit exists. Given f e 3^ (R), 
the (generalized) Fourier coefficients of f are defined as follows 

aA(f) :=M(f(t) e'^t) 

and the set 

cT(f) {AeR:QA(f)^0}, 

the spectrum of f, is at most countable. 

Define the equivalence relation ~ as f ~ g if aA(f — g) =0 for all A G R. As stated above, 
I • Isp is not a norm, as it can be zero on nonzero functions. It turns out that the quotient space 
■BP/ - is a Banach space (see [ ]). If f , g G 23^ (R) then the mean M(fg) is well-defined and 

M(f g) = M(fi gi) if f ~fi and g~gi. The space 23^/ is an Hilbert space when endowed with 
the scalar product (f, g)32 — M(fg) and we have the fundamental result due to Besicovitch. 

Proposition 3.2. The exponential functions t e'^* are an orthonormal Hilbert basis for 23^ (R). In 
different words, any f G 23-^ (R) can be represented by its generalized Fourier series 

(3.1) f(t)^ Y. 'lAme'^*. 

AGcr(f) 

fl«dX;^|aA(f)P < oo. 

Conversely, if one has a generalized series as above with square summable coefficients (as a generalized 
series), then there is a function f G 23^ (R) having the series as its own generalized Fourier series. 

Hence, we can identify a function by means of its generalized series as follows 

(3.2) 232(R):={fGL2„JR): ||f|||. := ^ |QA(f)P < +00} 

Ae(T(f) 

and the identification f (t) ~ ^Aea(f ) ) ^^^^ holds in the sense of convergence in 23^ (R). 

We turn to the pipe problem (2.2) in the almost periodic case. Since the problem is linear, it 
is reasonable to find a solution in terms of Fourier transform (or series). It is clear that, once 
the problem is solved in the Fourier space, we are given with coefficients a(f,) G C for £, G R 
and we are left with the problem of reconstructing the solution by inverse Fourier transforming. 
Since the convergence for classical Fourier series is robust in \? for coefficients, likewise in the 
context of almost periodic functions we consider the (correct) space 23^, for which the analogous 
of the Riesz-Fischer theorems holds true. 

In the following we shall also need spaces of the type Sobolev-Besicovitch, which are defined 
in the following way. 

Definition 3.3. Given a real s > 0, a function f G ^^(R) belongs to ■B^'2(R) if 

Wnls.:^ Y_ (l + |Ap) = |aA(f)P<+oo. 

Aecr(f) 

In particular, if f G 25^-^ (R), then the Fourier series for the (formal) derivative f of f is 
convergent and defines an element of 23^ (R). 

3.2. Construction of the "basic flow" in the Besicovitch setting. In this section we solve prob- 
lem (2.2) in the unknowns (w, tt), with the Besicovitch meaning. As a by-product of the method 
we obtain also a new proof of the existence of the basic flow in the periodic case. 
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Theorem 3.4. Given f e 'B^-'^[R) there are w e S^^R; H2(D)), with 9tW e ^^(R; L2(D]], and 

7t e S^fR) SMC/Z t/lflt 



9tW — vAw~7T 



and 



w(t,x) dx~f(t). 



Moreover, w and n are unique up to identification as almost periodic functions. Finally, there exists 
c > such that 

||Aw||32(L2) ^ c||f||cB2 + ;^||f'||32, 
||7t||s2 + ||9tW||252(L2(D)) ^ Cv||f||252 + c||f'||32. 

In order to prove the theorem, we restate the problem by spectral analysis in terms of Fourier 
transform with respect to the time variable t (with conjugate variable £,). Once Fourier trans- 
formed, problem (2.2) reads as follows: Find (w, n) such that 



(3.3) 



i£,w(£„x) -vAw{£„x) =7T(f,), 
w(^,x) =0, 
jDw{t,x)dx = f(^) 



xe D, £,£ R, 
X G ao, £, e R, 
f.e R. 



Clearly, the same result follows by a decomposition in Fourier series. The first equation yields 

w(£„x) = 7t(£,)W£,(x), where 

Wi-[x) := ((i£,Id-A/A)-il)(x) 
is defined to be the solution to the linear, stationary, and complex system 



(3.4) 

parametrized by £, e R. Set 



i£,Wi: -VAWj: = 1, 

W£(x) =0, 



X e D, 
X £ 9D, 



Wi^ix) dx, 



then by Fourier transforming the flux condition in (2.2) we get f (£,) = a£7T(f,) and in conclusion 
the solution to (2.2) (or, more precisely, to (3.3)), is given by 



(3.5) 



n[E,)^—n^) and w(£„x) = — fl^W^fx] 



The problem reduces to analyse the behavior of the two terms and with respect to £, e R. 
The main properties are summarised in the following lemma. 

Lemma 3.5. For every 1,^'R.it holds 

1. = V Jo |VW£(x)|2 dx - it Jo |W£(x)p dx, 

2. V Jo AWi:(x) dx = if,a£ - 1, 

3. v2 Jo |AW£(x)|2 dx + e Jd |W£,(x)P dx = 1. 

Proof. To prove the first property, take the complex conjugate of the equation satisfied by W^, 
multiply by Wf^ and integrate by parts, obtaining 



W£(x)dX: 



-itIWii- 



vW^AWfJ dx 



|VW£,(x)pdx-i£, 



W£,(x)r dx. 
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For the second property, just integrate the equation for on D. In order to prove the third, 
take again the complex conjugate of the equation for W^, but this time multiply by AW^^. Next, 
integrate by parts and use the first two properties to get 



i^a£ - 1 



which proves the equality. 



AWs^dx 



|AW£,pdx + iv£, 



iVWiP dx 



|AW£pdx + if,a£,-£.2 



|W£p dx, 



□ 



Next, we need to understand the growth/decay of and We^[x) with respect to £,, in order 
to show that the formal expression (3.5) defines a solution (in a suitable sense). 

Lemma 3.6. The map £, i-> is continuous on R with values in H^(D) n Hj(D). Moreover, as 
|£,| — > +00 zve have — ^ in (D) and also 



(3.6) 



lim £,|1VW£||2, =0, 

£,|— i-+oo 



and 



lim £,2||Wi:||[2 = 1. 



Finally, the map £, is continuous on R zvith values in R and 

lim a£, = — — . 

£,l->+oo £, 

Proof. Fix £,, £,o G R, with £, 7^ 0, and set V :— W^^^ — W^^. By symmetry we can assume that f, > 
and start with the case £,0 > 0. The new function V solves i£,V — vAV — i(f, — £,o)W£^q. Multiply 
by AV and integrate by parts to get 



it||VV||2, +^,||AV||2, =i(£-^o) 
The imaginary part of the above formula yields 



VWi:„ • VVdx. 



VW£„.VVdx) ^ |||VV||^L^ 



21, 



IVWi 



£,oIIl2 



and so ||VV||l2 < 



t.-£ol 
|£,l 



VW^f, 1 1 L2 . On the other hand the real part yields 



v||AV| 



2 

L2 



VWio • VVdx) s; 



s=:|£-tol||VW£j|L2||VV||L2 



1^1 ' 



and as f, — > f,o, continuity follows. 

In the case £,0 = the proof is slightly different, since one can prove directly taking the real 
part that 

^/||AV||22 ^ItlllVW^JI^^. 
Since V e Hq(D), this proves that as £, — > 0, then R V tends to zero in H^(D), but now with the 
order of !£,j^. 

Next, we consider the limit at 00. By the previous lemma we know that v||AW£^||l2 ^ 1, so 
there is a sequence W^^^ converging weakly in H^(D) to some Woo G H^(D). Indeed Woo — 0, 
since by the second property of Lemma 3.5 E,ai^ converges to a finite limit (and in particular 
implies (3.6)) and so 0, in particular Woo = 0. Moreover, from this it follows that the whole 
function W^ ^ weakly in H^(D). Finally, from the third property of Lemma 3.5 ||AW£^||l2 
converges strongly to 0. 

The statement on follows in the same way, by using again Lemma 3.5. □ 
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Remark 3.7. With a little more effort one can show that — > — |- point-wise as \E,\ +00, but 
we are not going to use this property in the sequel. 

Proof of Theorem 3.4. We use the identification (3.2) of an almost periodic function with its Fourier 
series and of the semi-norm with the sum of squares of Fourier coefficients (namely, Parse- 
val's identity). Consider f f(£,) e'^^, then we only need to prove suitable boimds for the 
quantities \ai^\, || AW£,||l2, and ||W£,||l2. 

Indeed, solving problem (3.3) for the Fourier components yields 

^(t, X) y e'^' and 7r(t) 2 V I^il e'^^ . 

In order to capture the dependence of constants from v, we observe that W^^ — ■^'W^/.y, where 
is the solution to (3.4) corresponding to y = 1. Set 

no[i)^\\W^\\l., ni(e) = ||VW£||^, UalO = ||AWt||^.. 

Indeed, 



since for |£,| < v. 



^ ci||f||32(jj,) + :^l|f'll32(jj,), 



=^ TTTTT ^ niax — -- = Ci, 



while for \E,\^ y, 

^'t12(|) ^ 1 ,2 Tl2(|) ^ 1 ^2 n2(ff) C2 



^ r?^ 77772 77-72 ^ 77^ sup — — = — ^ . 



v2ni(^)2 + t2^„(^)2 - v2 [(^)Vo(^)]' " v2 ,,^-1 (a2no{cT))' 
Similarly, 



2ni(^)2 + t2^o(f)2 



^C3v2|lf|||,+C4||f'|||., 



since for |£,| < v. 



V2^2^0(f) ^ V2 ^ 2 1 

E E — ^ r- ^ 'V max — -— = C3V , 



while for |£,| ^ v. 
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Finally, with similar computations. 



S^C5v2||f|||. +C4||f'|||2, 

where C5 = max|o-|^i nltW' "'"^^ quantities Ci, . . . , C5 are easily seen to be finite by the previous 
lemma. □ 

Remark 3.8. The computations in the proof of Theorem 3.4 provide an alternate proof to Theo- 
rem 1 in [ ], as well as to Corollary 2.9 (once (generalized) Fourier transform are replaced by 
Fourier series). Moreover, the following estimates hold. 



|Aw||[, dt^c||f||2.(^,+-||f'||2.(j,,, 



9tw||L2 dt - 



|7r(t)pdt«;cv2||f||2.,^,+c||f'||2.,^,, 

R 

^ |2 



sup ||Vw(t)||2, ^ C(v||f||2, + -||f' 

teR Y 
Indeed, in the proof above we have shown that 

MM!i^cmax{l,v-'|!,|} 

< cmax{|£,|,v}. 



Hence, by Parseval's identity, we get 



IIAw||l2 



^ M)?\\^yv^\\h^m\i + ^\\rf 



— — ^|HC,J| ||^VVi:||L2 ^ L||1||L2 -Jill ||L2, 

R l^il V 



JR 

and the other inequalities are obtained similarly. Finally, the inequality for Vw follows by inte- 
gration by parts and the identity \ \ Vw| | ^2 = —2 Wt Aw dx. 

3.3. On the meaning of the solution. We need to spend a few words about the notion of so- 
lution we constructed. We observe that a given f e ®^'^(R) is clearly identified in the sense 
of ®^'^(R), hence by means of its generalized Fourier series. This implies, for instance that 
if w is a solution in the sense of Besicovitch spaces, then w + w is also a solution, for any 
w e L^(R.; H-^) n H^(R; L^). This poses some restrictions to the interpretation of the result. One 
would like to have some embedding in the space of continuous functions in order to have a 
more precise identification of the solution. A larger spaces in which we are able to solve the 
equation is balanced by a weaker notion of solution. 

In general one cannot expect the validity of the usual Sobolev embeddings in 2?'^'P(R) as is 
explained for instance in Pankov [ ] and especially the identification with UAP(R) functions 
is not a trivial fact. Classical counterexamples can be found in the references cited, while the 
following general embedding result is proved for instance in [ ] . 

Proposition 3.9. Let E c R ^^e countable and assume there fs |3 > such that the generalized sum 
satisfies 

I J < +00 /or y > (3 
f^J^ |=+oo /ory < (3. 
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J/|3 < 2s, then for every f e 'B''^'^[R) such that (j{f) c E, lue have f e C^'"(R) n UAP(R)/or all 
a e [0, s — r — |3/2], ivhere r — [s — !■] (ivith corresponding inequality for the norms). 

Remark 3.10. To simplify the notation from now on we denote by f ^ the (generalized) Fourier 
coefficient f(f,) that is (more precisely) written as a£,(f) in (3.1). 

To understand this result, let us observe that if G ^^{C], then the series e'^* con- 

verges uniformly and can be identified with a continuous almost periodic function f ^ C£^ e'^*. 

Moreover, for classical Fourier series, i. e. ff(f ) = E C Z, the (3-condition is satisfied for P = 1 
and this shows that if (cj)jgz/ (jCj)jgz G f^(C), then 

jez iez iez jez\o 

This is the guideline to understand the result for 23'''^ (R), since one has - roughly speaking - to 
show an inequality similar to 



£,ecr(f) f,ecr(f) £,ea-(f) £e(T(f) 

For instance Proposition 3.9 implies the following result. 

Corollary 3.11. Letbegiveni e ■B^-'^(R), zvith i-^^^^s^e'^^ such that 

Then, there are w e 'B^[R;H^{D)) n UAP(R; HjlD)], with dtw e 'B^[R;L^[D]), and n e 'B^{R] 
such that (2.2) is satisfied in the sense of Besicovitch. 

This makes also possible to consider the flux as 

f(t) =fi(t)+f2(t), 

with fi e 'B^'^[R) and f2 G H^(R), so that ||f2||3i,2 — 0. One can construct the solutions 
wi e 23^'^(R; L^(D)) and W2 e W^'^(R; L^(D)) corresponding to fi and f2 respectively, and add 
together. 

This is not completely satisfactory, since we still do not have a precise identification on the 
pressure. To this end one would like to have a solution in the classical UAP(R) space for exam- 
ple. This can be achieved by assuming stronger conditions on f (rather than on its spectrum), 
as shown by the following result. 

Proposition 3.12. Let f be given, with f ~ f £, e'^^, such that 

^(i + itnif£i < +00. 

Then, there exists a unique solution (w, tt) to (2.2) such that 

(3.8) ^ ^ ^ 

IK II 
where w^^ and are the (generalized) Fourier coefficients ofw and n, respectively. In particular, w e 
UAP(R; H2), 9tw e UAP(R; H) and n e UAP(R). 
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Proof. First we notice that since |ff,| < +00, then 

showing that f G 23^ (R). The same argument shows also that f e 2?^(R), and so Theorem 3.4 
ensures the existence of a unique solution. Since wi^ — a^^^W^ff^ and tt^ = a^^^f^,, the es- 
timates (3.8) follow immediately from (3.7). In order to show that w and n are Bohr-aLmost 
periodic, we consider a truncation 

where (crn(f))TieN is an increasing sequence of finite subset of cr(f) such that IJ^ cr^(f) — cr(f). 
For each n e N we can consider (2.2) with flux given by the trigonometric pol5momial fn and 
the estimates (3.8) imply uniform convergence of the corresponding solutions (wn, Tin) towards 
UAP functions with the requested properties. □ 

3.4. The nonlinear case. In this last section we consider the non-linear problem. Assume pre- 
liminarily (we shall assume stronger assumptions on f later) that f e 23^'^ (R) and denote by Wi, 
W2 the basic flows in the two pipes Oi, O2 respectively, provided by Theorem 3.4. Let Vq be the 
flow defined as in (2.21), it is clear that Vq is also almost periodic and keeps the same regularity 
properties of wi and W2, namely 

Vo e S2(R; h2(0) n HliO)) n 'B^'^{R; L'^iO)), 

as well as the flow w defined in (2.22). Indeed, both flows are obtained by applying only linear 
operators in the space variable to Wi and W2. 

Consider the full nonlinear Leray's problem in the (Besicovitch) almost periodic setting, 
namely to find a solution (u, p) to the problem 



j 9tu - vAu + (u • V) u -h Vp ~ 0, 

||U- Wi||B2(H,.Hi(Oi)) s% Co, i = l,2, 



(3.9) 

such that 
and this implies that 

(3.10) lim ||u-Wi|l232(R.Hi/2(Di)) = 0- 

If w is the flow defined in (2.22), consider the solution u = U + w as a perturbation of w. 
Consequently, 



(3.11) 



' 9tU - vAU + (U • V) U + (U • V) w + (w • V) U + VP ~ F, 
V-U~0, 



where F is defined in (2.24). The main theorem of the section is the following. 

Theorem 3.13. Assume that the flux f satisfies 

(3.12) 4), :=^(l + |t|)|fd <oo. 

Then, there exists Vq > 0, with Vq — vo(f, O), such that for every v ^ Vq problem (3.9) admits a 
solution 

ue S^IR; Hj(0))rTB5'2(R; \}[0)), 

which satisfies (3.10). 
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The rest of the section is devoted to the proof of this result. 



3.4. 1 . Spectrum and module. Before turning to the analysis of problem (3.11), we recall that, since 
f e (R), its generalized Fourier series is well-defined and its spectrum (j(f ) is the set of modes 
£, e R corresponding to non-zero coefficients f £^ in the Fourier expansion of f . 

Since f is real, it follows that — f_£^ and so the spectrum is symmetric, namely — cr(f ) — 

(T(f). 

Definition 3.14. The set |a.(f ) is the Z-moduIe of the spectrum of f, namely the smallest subset of 
R which contains (T(f ] and is closed for the sum (that is, if ti e M.(f ], then a £, + b ti e M.(f ), for 
all a, b e Z). 

It is clear that ^(f) is also symmetric and, since cr(f) is at most coimtable, |J.(f) is at most 
coimtable too. Moreover, it is easy to see that the spectra of w and F, by linearity, are contained 
in the spectrum of f . Indeed, by construction, the terms Vo and w, defined in (2.21) and in (2.22), 
respectively have spectrum contained in cr(f ). 

In the following, with the purpose of approximations, we shall need to consider finite dimen- 
sional truncations. To this aim, we fix an increasing sequence (|j.n (f))NGN of subsets of [a.(f) con- 
verging to n(f), that is M-N (f) C M.N+i(f) and = Un (f), and such that ^N (f) = -M-n (f)- 



3.4.2. Reduction to a system in Fourier variables. A remarkable feature of the nonlrnearity we are 
going to analyse is that if ui, U2 are Besicovitch almost-periodic, then (ui • V) U2 is also in 
the same class. This result on product of almost periodic functions is not true in general, but 
as we will see in our case it holds since the spectrum of the nonlrnearity is contained in the 
module generated by cr(ui ) and cr(u2 ) . Having this in mind, we recast problem (3. 1 1) in Fourier 
variables, 

(3.13) itU£-vAU£+ Y_ [(U^-V)Ue + (U^-V)we + (w^-V)Ue] +VPt = ft, 

Ti+e=£, 

for £, e H{f ), with V • U^^ = 0, and the sum in the formula above is extended over all r], 9 G H(f ). 

We shall use the following strategy to prove Theorem 3.13. We linearise the nonlrnearity (by 
introducing an auxiliary field U) and solve the new linearised problem (in two steps, first for 
a finite number of modes, then for all modes). The assumption on the viscosity allows to have 
a uniquely defined map that gives a solution to the linearised problem for each field U. The 
same assumption ensures that this map is a contraction and its fixed point is the solution to 
problem (3.9). 



3.4.3. Preliminary tools. We prove two preliminary tools for the analysis of the problem. 

We first consider the fields w and F defined respectively as in (2.22) and (2.24) and prove the 
following estimates in terms of (})*. 
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Lemma 3.15. Let be given f e !B^'^(R), assume that (3.12) holds. Then, there is c> (independent of 
v) such that 

||Aw||s2(R.L2) < c(l + :^)||f||3i.2(R), 
I|7T||s2(r) + ||9tW||T32(R;L2), ^ C(l +A')||f||a3i,2(R) 

^IIAw^IIl^ <c(l + i)4)^, 

||F||32(R.L2) S; c((l+A/) + (l + i)2ct).)||f||3i.^(R), 
^||f£||L2 <c((l+v) + (l + i)24).)4).. 

Proof. The inequalities for w, tt are a straightforward consequence of Theorem 3.4, Proposi- 
tion 3.12 and the definition (2.22). To prove the inequalities for F, we only need to consider the 
term (w • V) w (the estimate of the other terms in F follow from the estimates for w and tt). For 
f, e M-(f ), writing the explicit expression for [[w • V) w]^ = ^Ti+e=£,('^Ti ' ^) obtain that 

||[(w V)w]^||l2 < c Y_ l|Vw^||L2||Awe||L2, 



hence 



(w • V)w|||2(R.L2) s$ Y_ l|VwTi||L2||Awe||L2 



£, ri+e=£, 

N 2 

£.IIl2 ) ||^W||252(I{.L2 



^ ( > liVwdlr^ 1 ||Aw"2 



and also 



^||[(w- V)w]£,||l2 Y_ ||Vw^||L2||Awe||L2 

£. £, Ti+e=£, 

which complete the proof. □ 

The second result is the extension to our setting of the usual cancellation property of the 
nonlinear convective term, when energy estimates are derived. 

Lemma 3.16. Let be given X = ^ X^, e'^* andX ^ ^'^^ ^o*^ belonging to 'B^[R; Hj), with Y 
divergence- free. Assume that the spectra ofX and Y are contained in [i[f) and fix an integer N ^ 1. 
Then 



^(N) y- (N) 



X£ • (Y^ • V] Xe dx = 0, 



where the superscript on the sum above means that the sum is extended only over modes in [J-n (f )• 
Moreover, the same holds true for N — oo if at least one between X and Y is in 2?*(R; Hq), with 

a3*(R;Hi):={feLi„,(R): ||f||s.(R,Hj) := L (1 + l£,l)l|fdlHi< +oo}. 
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Proof. Let us denote by @ the sum in the statement of the lemma, then by a change of summa- 
tion index 



^ (N ) f ^ (N) 

X£.(Y^.V)Xedx = -^ 

£,+Ti+e=o ° £,+Ti+e=o 



Xe • (Y^ • V) X£ dx, = ~@, 

o 



since Xf^ = X_£, and the claim is true. To show that the same holds for N — oo, it is sufficient to 
prove that the following sum extended is bounded uniformly in N . Indeed, if Z is another field. 



VYT,||L2||V2e||L2 



£, ri+e=f. 



• [Y^ • VJ Ze dx ^ } 

o 



^ ||VX||;B2(h1)||VY||32(h1)||Z||;B*(HJ)i 

by Young's inequality for convolutions. □ 

We observe that in the cancellation property above it is fundamental that we deal with the 
complex conjugate of X and with the fact that the spectrum is symmetric, since flux and solution 
are both real-valued. 

The proof of Theorem 3.13 is split into three preliminary steps. 

3.4.4. First step: existence for the finite modes approximation. Given U e 23^ (R; Hj) with V • U = 
and an integer N ^ 1, we seek for a solution to the following problem, 

(3.14) iEME, - vAUt + [(U^ • V)Ue + (U^ • V)we + (w^ • V)Ue] + VP^ = f^, 

Ti+e=f. 

for £, e [O-N (f ), where again the superscript on the above sum means that the sum is extended 
only over modes in (f)- 

Proposition 3.17. Let f e B^-^tR) and U e ^^'^(R; Hj) and set 

4)^ ^'^'(1 + If.l)!f£! and ^t^' 

Then, there is at least one solution U''^' to problem (3.14). Moreover, there are non-negative Ci(-, . . . ) 
and C2(-, ....... ) increasing functions of their arguments and (depending only on the domain O) such 

that 

(3.15) ||U(n'|1^..,j^^^,j + ||U(n'||b.(R;hi) <c(vJ|fbi,2,K,)xl;^+c(v,4,^||f|h...(K^ 
If additionally 

v>ij;r +c(l + i)4>^ 
where c is the constant of Lemma 3.15, then the solution is unique and 

(v-^j^r-c(i+i)ct,r)(x'"iivu<^'ii,; 

(3.16) ^ £ 

^ c(l + y)4>^ + c[l + i)2(4).^ f + c(l + i)4>rtJ^r • 

Proof. The proof can be carried on with the standard technique of Fujita (cf. Theorem 1.4 of [27, 
Ch. 2]) for the case of existence of solutions for the steady Navier-Stokes equations in un- 
bounded domains (but we have the additional advantage of the Poincare inequality (2.20)). We 
use Galerkin approximations — ^.^=1 '^^^ic (not necessarily made with eigenf unctions) and 
we consider the projection of (3.14) on the finite dimensional Galerkin space as a problem on 
jj^2nN ^|.]^g j.g^j gj^j imaginary parts of each u'^ count as two variables) with the scalar product 
induced by the one of Hj (0)®2N . 
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We show existence of a solution of the finite dimensional problem by means of Lemma 1.4 
of [27, Chapter 2] (which in turns is a consequence of Brouwer's fixed point theorem). Let P(ll^) 
be given component-wise by the projection of (3.14) so that if P(U^) — then is the solution 
to the Galerkin projected problem. It is sufficient to show that P(U"^) • Un > on |U^| = Co for 
some Co > 0, where product and norm are those we have given on R^"^"^ . This is immediate by 
Lemma 3.16 since 

P(U-) • U- > v|U-p - YT^ |U£|||VU^||L2||Vwe||L^ -c||F||3^(L^)|U-| 

£ r|+e=£. 

;^vlmp-^l,N||Vw|ls2(L^)lU-|-c||F||,£.(L2)|U-|, 
which is strictly positive if we choose v such that 

Co > ^(lj)I'||Vw||s2(L2i +c||F||s2(L2)). 

Passing to the limit in the Galerkin approximation is standard (the non-linearity contains a finite 
sum) and follows from uniform bounds (in n) on U*^ which are similar to (3.15) and whose proof 
is formally similar. Hence, we prove (3.15) directly. For each £, e |a.N (f) multiply (3.14) by U^^, 
integrate by parts on O, sum over £, G m-n (f )/ and use Lemma 3.16 to get 

X— (N) n X— (N) 



(N) ^ (N) ~ V— 

£, Ti+e=£. ■'o £, 
Flence by Young's inequality for convolutions and taking the real part we get 

(N) ^ (N) ~ 
V||VU|||.,L2) ^ 2_ 2_ ||VU£||L2||VU^||LH|Vwe||L2 +C||F||3(L2)||VU||32(L2) 

£, Ti+e=£, 

l|;^||Vw||32(L2)||VU||s2(L2) +C||F||3(L2)||VU||32(L2), 

and so inequality (3.15) follows from Lemma 3.15. 

To prove (3.16), multiply (3.14) by U^, integrate by parts on O and divide by the non-zero 

||VU£||l2 to get 

v||VU£||l2 c||f£|| + ^'^'(||VU^||L2||VUe||L2 + ||VUJ|L2||Vwe||L2 + ||Vw^||L2||VUe||L2). 

ri+e=£. 

Inequality (3.16) follows by summing in using Young's convolution inequality and Lemma 3.15. 

Finally, if U^*^', \^^ ^ are two solutions corresponding to the same data, let D'*^' := uj'^' — 
U^^'andQ(^':= pj^' - P^^', then 

itD^^' - vAD^^' + X""' (Un • V) D^^' + (w, • V) d[,^') + VQ^^' = 0, 

ri+e=£ 



(N) 



O 



and hence taking the scalar product with D ' 

V||VD'^'|||2,K;L2, ^ ||VW,|k2||VDf^'|k2||VD[,^>|k2 

£, ri+e=£, 

^(X'^'||Vw,||l2)||VD(^'|||2,k;l2„ 
which, by the assumption on v, implies that D = 0. □ 
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3.4.5. Second step: existence of a limit as U oo. Letf e 23^'^(R)nS*(R) and assume additionally 
that the quantity 

is finite. This implies in particular, as in Lemma 3.12, that f and f ' have representatives which 
are Bohr-almost periodic. Let U e 53-^ (R; Hq), assume that the quantity ||VU£||l2 is also 
finite and consider the problem 

(3.17) itU£-vAUt+ [(U^-V)Ue + (U^-V)we + (w^-V)Ue] +VP£=ft, 
for t e ^(f). 

Proposition 3.18. There exists Vq > 0, with vq = Vo(f, O], such that for every v ^ vq there is \|;^, > 
such that 

V >tl)^ + c(l + i)4)^, 

and if 

then there is a unique solutionUe ^^(R; Hj(0]] n 2^/2,2^5^. l2(0)) to problem (3.17). 
Moreover, 

ll^llBi'^(R;L2) + I1U|1b2(R;H1) =^ c(v, || f|| 231^2 (r, + cl^/, , ||f || SI.2 (r) ) , 

(3.18) _^||vUt||L2 

Proof. Let (U''^')n^i be the sequence of solutions to (3.14) provided by Proposition 3.17. By (3.15) 
it follows that (U(^')n^i is bounded in S^/^'^fR; L^) and in ^^(R; Hj), hence there is a sub- 
sequence weakly convergent to a limit point U G "B^^^'^IR; L^) n ^^^(R; Hj). Since (3.17), weak 
convergence is enough to pass to the limit in the equation. Uniqueness follows as in Proposi- 
tion 3.17, using the bound on the viscosity. 

We only have to identify Vq and ij;*. From (3.16) it follows that 

{y-^\>,- c(l + i)4).) {Y_ II VUdk^) < c(l + y]4>. + c{l + + c(l + ^]<p.^\>., 

£, 

SO everything boils down to show that for v large enough there is x[>* such that 

C(l + v)cl), + c(l + ^)2(4,^p + c(l + i)4).^. ^ , 

1 ^ w*, 

that is 

^\>l-{y- 2c(l + i)ct).)tl^* + c(l + v)4). + cHl + i)'*' ^ 0. 

It is elementary to verify that the above polynomial has two positive solutions for "v large 
enough. □ 

Remark 3.19. Clearly, without the assumption on the size of v in the previous proposition, one 
can still show existence of at least one solution to (3.17). The size condition on y is necessary 
only for proving uniqueness. 
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3.4.6. Third step: the fixed point argument. Under the assumptions of Proposition 3.18 we have a 
well defined map U U, where U is the solution to problem (3.17). Denote the map by J, then 
it is clear that any fixed point of J is a solution to (3.13) and hence to (3.11). 

Proof of Theorem 3.13. Fix y ^ Vq, where Vq is given in Proposition 3.18. We prove that the 
map J is a contraction on the set X of all U e 'B'^[R; Hj(0]) n 'B^/'^^'^[R; L'^[0)) that verify the 
bounds (3.18). 

The fact that 3 maps X into X clearly follows from Proposition 3.18, so we only need to prove 
that 3 is a contraction. This is obtained as in the proof of uniqueness of Proposition (3.17). 
Indeed, if E := Ui - Us and D := J(Ui) - 3{U2), then 

itE£-vAE£+ Y_ (E^ -VJUi +(U2 .V)De + {E^ -Vlwe+VQi^O, 

Ti+e=£ 

with a suitable Q. By multiplying by D^^, integrating by parts, and summing over £, we get 

(v-tJ;,)||VE||s2(L2) ^ (^, + c(l + ^)4),)||VE||32(L2) 

that is ||VE||32(L2) ^ Ko|| VE||32(l2), with 

^. + c(i + i)4)^ 

Likewise we also have i aj^^sj ^ "^1^oII^E||32(l2). Finally, by multiplying by D£_, dividing 
by 1 1 VD £^ 1 1 L2 and summing over f, we get 

f. £, 

In conclusion, if v is large enough it follows that Kq < 1 and vKq < 1 and the map 3 is a 
contraction. 

In order to conclude the proof, we need to show that if (U^j^^g^ff) is solution to (3.13), then 
the Besicovitch almost periodic vector field U having Fourier coefficients (U^j^^g^ff) is a weak 
solution to (3.11), namely that for every divergence-free cp e C^(0; R^) and every £, e R, 

M ( (9tU, (p> - v(U, Acp) - (U, ((U + w) • V) cp) - (w, (U • V) cp)) e'^^l = 0, 



which is an easy consequence due to the bounds (3.18). □ 

3.4.7. Final considerations. Apparently the assumption (3.12) seems to be essential for the proof 
in the Besicovitch setting to work. The technical problem is essentially related to the term 

£ Ti+e=£, 

which is of order three, although all bounds on U, V, W are of order two, if one works in the 
framework of 25^ spaces. Young's convolution inequalities tell us that in general there is no 
possibility to bound the above term under these assumptions. In terms of the time variable, we 
are trying to bound the Navier-Stokes nonlinearity over the whole R. 

Another possibility would be to use the other a-priori estimate, namely the bound in 23 ^^^'^ (R; ), 
which plays no role in the proof of Theorem 3.13, using for instance the results in Section 3.3. 
This possibility is rided out by the non-linear term. In fact, in the standard case of Leray-Hopf weak 
solutions one has a better knowledge of the time derivative and this can be used for instance 
with the Aubin-Lions compactness lemma to handle the non linear term. 

Indeed the non-linearity reads in Fourier variables (in time) as a convolution and, whatever 
is the spectrum cr(f) of the flux, the spectrum of the solution to the non-linear problem will 
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have the Z-module |a.(f ) as its spectrum. In different words, the non-linearity creates a full set of 
harmonic resonances in the time frequency. The structure of Z-modules in R shows that the only 
possibility to use a bound on the derivatives (while obtaining a useful information for all times 
t G R) is the periodic case, previously studied in [ ] . Indeed it is easy to verify the following 
result. 

Proposition 3.20. Let G cRbe a Z-module. Then, either G = k Z/or some k g R or G zs dense in R. 
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